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Abstract 

We present observations of the closest region of high 
mass star formation, Orion BN/KL, performed at both 
low resolution mode (grating mode) and high resolution 
(Fabry-Perot mode) with the Long Wavelength Spectrom- 
eter (LWS) on board the Infrared Space Observatory (ISO). 
We detected the CO rotational transitions from Jup— 15 
to Jup = 49. A LVG analysis of the line fluxes allows to 
distinguish three main physical components with different 
temperatures, densities and column densities. Our analy- 
sis yields to conclude than J up <15 arise from the pho- 
todissociation region (PDR), the emission between J„p= 
20 and 30 arise from the high velocity outflow (plateau), 
and Jup > 32 arises from a hot and dense gas component. 
The latter exhibits broadened lines for the levels Jup > 32 
and is thought to be due to shocked gas in a high veloc- 
ity outflow. Future observations with HIFI, onboard the 
Far Infrared Space Telescope (FIRST) will allow the spec- 
tral separation of the PDR and the plateau components, 
unresolved with ISO, and characterise more precisely the 
Orion BN/KL star forming region. 



1. Introduction 

The Orion Molecular cloud, at a distance of 450 pc, is the 
closest region of high mass star formation. Its proximity 
and its large infrared luminosity allowed to perform plenty 
. of observations in the past years, yielding the discovery of 
' the first proto-stars candidates. Molecular emission mostly 
comes from OMCl, containing several condensation as the 
KL nebulae. The KL nebulae is composed by many in- 
frared clusters (e.g. BN or IrC2), containing massives stars 
at early evolutionary states. Millimeter, sub-millimeter 
and infrared spectroscopy fsee lGenzel fc Stutzki 19891 for 
a review) have shown copious molecular emission arising 
from physically distinct regions: the ridge, the compact 
ridge, the hot core, the PDR region surrounding the qui- 
escent gas and a high velocity (Aw = 18 km/s) bipolar out- 
flow originating from IrC2. Outside the bipolar outflow is 
a region of very hot (1000 to 2000 K) shock exited gas. Nu- 
merous observations have allowed to charaterise this com- 
ponents in term of temperature, density and column den- 



Spectrometer (hereafter LWS: |Clegg et al. 2000|) on board 
of the Infrared Space Observatorv fISO: IKessler et al. 1996|l . 
Theses observations revealed the presence of three gas 
components they identified as the ridge, the high veloc- 
ity outflow and a hot and dense gas component detected 
at Jup > 30 due to the shocked gas in the high veloc- 
ity outflow. Notheless uncertainties remains on the data 
calibration, and high J grating data are probably contam- 
inated by adjacent lines. We here present a large bandwith 
survey of the far infrared CO lines of Orion BN/KL, per- 
formed with ISO-LWS in Fabry-Perot mode. Particular 
emphasis is given to the calibration of the FP data. Based 
on theses new calibrated data we interprete the FIR lines 
by means of a LVG code to derive the temperature, den- 
sity and column density of the several gas component. 

2. Observations and results 

We performed a spectral survey of the Orion BN/KL us- 
ing ISO-LWS both in grating mode and FP mode. The 80" 
beam was centred on Orion BN/KL (a2ooo = 5''35™14.2^ 
(52000 = — 5°22'33.6"). The grating spectral survey was 
done using LWS in with the LOl AOT. It was calibrated 
using Uranus, and the absolute accuracy is estimated to 
be better than 30% ( |Swinyard et al. 1998| ). Theses obser- 
vations required the use of the bright source data reduc- 
tion package IjLeeks et al. 1999j) as the brightness of the 
source saturated the detectors LW2 to LW4. The grating 
spectrum was only used to calibrate the FP observations. 
The FP spectral survey was performed with the L03 AOT. 
Theses observations, covering the wavelength range 43 to 
162 fim are the first performed with a single instrument 
in space at the same time. The continuum level of the 
observation was calibrated against the grating spectrum, 
after taking into account the dark, straylight and FP or- 
der sorting. The calibration uncertainties, estimated from 
RMS of measured fluxes of each transition, is between 30% 
and 40% for the weaker lines. We detected ^^CO rota- 
tional lines between Jup = 15 and 49 in FP mode. Higher 
transitions fluxes are under 10^^^ erg.s^^.cm^^ and are 
not detected. The fluxes uncertainties, estimated from a 
measurements of each transitions, are about 30% for tran- 



sitions between Jt, 



15 and 40 and 40% between 40 and 



sity. Recently Sempere et al. (2000) reported of FIR CO 
observations towards BN/KL with the Long Wavelength 



45. Lines between Jup = 45 and 49 must be considered 
as upper limits. Higher transitions have a low RMS yield- 
ing uncertainties on width determination but the broaden- 
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Figure 1. Observed profiles of selected CO rotational 



transitions whit LWS in FP mode, 
shows an important broadening 



Juv > 32 transitions 



Figure 2. Observed CO rotational transitions by LWS 
in FP mode (crosses with error bars), and three compo- 
nents model of Orion BN/KL. Uncertainties on LWS FP 
lines fluxes are given at 1 a. We also plotted the values 
from anteriors observations. The differences between this 
mesurements and LWS mesurement are probably due to 
the differents instruments beams. PDR emission is rep- 
resented in blue, plateau emission in orange and shocked 
gas emission in red. The total emission is represented in 
red. Predominance of the shocked gas component from J^p 
= 32 allows to explain the observed broadening at these 
transitions. 



Jup = 16 and the "broad shaped" emission between J„p 
= 20 and 30. In addition the observed broadening at 3 up 
> 32 impHes a third physical component. 



3.1. Low J EMISSION. Ridge, compact ridge and PDR 



ing observed is consistent with higher transitions measure- 
ment. This is clearly indicate that J„p > 32 lines emission 
comes from another gas component than lower transitions, 
necessary hotter and denser. We also detected J„p — 18 Because of the accurate calibration we used, the emission 
and 24 ^^CO lines. We found a flux of respectively 3.7 ± 
1.1 and 4.7 ± 1.4 erg s"^ cm"^. The 



ratio is 27 for Jt, 



18 and 12 for J 



fluxes 
24, indicating 



that theses lines are optically thick. 



up 



3. Discussion 

We analysed the ^^CO lines fluxes by means of an LVG 
model developed by Ceccarelli et al. This model, which 
compute in a self consistent way the opacities of lines, 
has four free parameters : the ^^CO column density, the 
H2 density, the gas temperature and the linewidth. The 
results of our computations are shown on fig. |21 together 

with our observations and previous ones IjSchultz et al. 19921 / 20 ^^CO lines ratios. This ratio, both with the photon 
[Graf et al. 1990; 'Schmid-Burg k et al. 1 989','Ho we et al. 1993fc scape probability were used to constrain the temperature 

and density of the gas. The fig. 13 shows theoretical pho- 



between J^p = 15 and 20 we observed is quite different 
than previous ISO observations. The fluxes we measure are 
30 % higher than those measured by |Sempere et al. 2000| 
Such fluxes can't be explain by the ridge emission as they 
previously said. The temperature, density and column den- 
sity they adopted fail to reproduce the intensity of the 
emission we observed. In order to characterise the gas 
emitting at these transitions, we used the ^^CO to ^^CO 
line ratio at 3 up = 18 to calculate the escape photon escape 
probability at this transition. Assuming that the relative 
abundance of ^^CO with respect to ^^CO is 60 and that 
the '^•^CO is optically thin, it gives a photon escape prob- 
ability of 0.5 at 3up = 18. We also calculated J„„ = 15 



up 



[Gcnzel et al. 19 88). The differences between these previ 
ous measurements and our measurements are probably 
due to the different instruments beam. The first thing to 
note is that a single gas component can not explain the 

3up = 15 to 44 observed emission. At least two compo- density of 10^^ cm^'^. We adopted a column density of 10 
nents are needed to explain the emission peak observed at cm~^ and a line width of 10 km/s. This parameters gives 



ton escape probability and lines ratios computed with our 
LVG codes, both with photon escape probability and ratio 
we observed. We found a lower limit for the CO column 



■,18 




Figures. Theoretical Jup =15/20 ^^CO lines ratios (in 
blue) and photon escape probability at J^p = 18 (in black) 
predicted by our LVG model. The observed values are rep- 
resented by continuous lines and the errors bar by dashed 
lines lines. We assumed N(CO) = 10^^ cm^^ and Av = 
10 km/s. This parameters gives a lower limit of 250 K for 
the temperature and 1(P cm^^ for the density. The best fit 
model is obtained for n(H2) = l(fi cm~^ and T = ^00 K. 

a a temperature lower limit of 200 K. The best fit model 
is obtained for T = 350 K and n(H2) = fO^ cm'^, but ob- 
served emission can also be explain by lower temperature 
and higher density, or inversely. The lack of data between 
J«p = 7 and 15 not allow to constrain more precisely the 
physical parameters. The parameters we adopted require a 
beam dilution factor of 1, which implies an extended emis- 
sion. This, both with temperature limit of 200 K we ob- 
tained, indicates that the emission between 15 < J„p < 20 
can not only arise from the ridge, but may arise from a the 
PDR region, in agreement with lHowe et al. 19931 Nothe- 
less, BN/KL is a complex region where a lot of gas compo- 
nent whit nearby physical characteristics are present. This 
components are not resolved with LWS. Even if a single 
gas component can explain the J„p = 15 to 20 emission 
we observed, a part of emission may arise from the ridge. 
Nonetheless, because of it low temperature, we estimated 
the low J emission contribution of ridge of 10 % of the 
PDR emission. The hot core, due to its small angular size 
towards the PDR, also certainly few contribute to the 15 
< J«p < 20 emission. The density and column density are 
in agreement with previous works (|Howe et al. 1993jl . but 
the temperature we derive is significantly higher. 

3.2. J„p = 20 TO 30 EMISSION. Plateau 

Observed emission between J^p = 20 and 30 show a broad 
shaped emission necessary arising from an hotter and denser 
component than the PDR. Along the same lines, we used 
the -^•^CO to ^^CO fluxes at J„p = 24 to constrain the 
physical parameters of the gas emitting at theses transi- 
tions. This ratio give a photon escape probability of 0.2. 



Figure 4. Theoretical J^p = 26/20 ^^CO lines ratios (in 
blue) and photon escape probability at Jup = 2^ (in black) 
predicted by our LVG model. We assumed N( CO) = 5.10^^ 
cm~^ and Av = 30 km/s . It gives a lower limit of 5.1(fi 
cm^^ for the density and 300 K for the temperature 



We used the ^-^CO line ratio 20/26 is 1.4 and adopted a 
line width of Av = 30 km/s. 

We found a lower limit for the column density of 10^® 
cm~^, but higher density are also giving solution. We 
adopted N(CO) = 5.10^^ cm^^. This column density gives 
a lower limit of 5.10^ cm^^ for the density and 300 K for 
the temperature (see fig. EJ. The parameters which are 
best fitting our observations are T = 350 K and a n(H2) 
= \{f cm~^. Theses requires a beam filling factor of 7 % 
to fit the iup = 20 to 30 lines fluxes. The angular size 
we derive from this filling factor is 20 " IHowe et al. 19931 
The high velocity component angular size deduced from 
CO 7-6 line emission is 40 ", but could be smaller at higher 
transitions. The angular size we determined yields to say 
that the emission at theses transitions arises from this 
gas component. The physical parameters we derived agree 
with previous determinations (,Schultz et al. 199'5)l . 

3.3. High J emission. Shocked gas 

Although high J emission could be explain by higher plateau 
temperatures and densities, the observed broadening of 
Jtip > 32 lines clearly shows that emission at this transi- 
tions arises from another hot and dense gas component. 
High temperatures and densities that requires this emis- 
sion suggest that the emission arise from the gas shocked 
by the high velocity outflow, in agreement with |Sempere et al. 2000| 
interpretation. We modelled this emission by a gas tem- 
perature of 1500 K, n(H2) = 4.10^ cm^^, and Av = 50 
km/s. Theses requires a beam averaged column density 
of 10^^ cm"^. Assuming a angular size of 40 ", it im- 
plies a column density of 10^® cm~^. This values are in 
good agreement with |Sempere et al. 2000| At J^p > 32, 
the shocked gas emission become more important than the 
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Figured. Jup — 18, 24- and 33 CO rotational transi- 
tions, superposed on PDR emission (blue), plateau (or- 
ange), shocked gas (red), and total (green) predicted by the 
LVG model. Theoretical line profiles have been obtained in 
convolving gaussians profiles of respectively 10, 30 and 50 
km/s width and of intensity deduced by the LVG model, by 
the PSF of the instrument. Predominance of shocked gas 
emission at Jup — 32 can explain the broadening of lines 
observed at theses transitions. 



plateau emission, explaining the line broadening observed 
at theses transitions. 

3.4. Agreement with observed lines profiles 



Table 1. Temperatures, densities, column densities and mi- 
cro turbulent velocities of ridge, PDR, plateau and the 
shocked gas determinated from our LVG model 





PDR 


Plateau 


Schocked gas 


n(H2) cm-^ 


10" 


2.10" 


4.10" 


T (K) 


400 


350 


1500 


N(CO) (cm-2) 


10i« 


S.IO^^ 


10^^ 


/S.V (km/s) 


10 


30 


50 


Filling factor (%) 


100 


7 


25 



4. Conclusions 

The Orion BN/KL observations with LWS allowed to de- 
tect J„p = 15 to 49 CO rotational transitions in Fabry- 
Perot mode. The modelling of observed fluxes by a LVG 
showed that molecular emission can be explained by three 
gas components that we characterised in term of temper- 
ature, density and column density (see tab.^l. We found 
that low J emission can not only arise from the ridge, too 
cold, but may arise from the PDR region. The emission 
between iup — 20 and 30 certainly arise from the high 
velocity outflow. Finally, observed J„p > 32 lines broad- 
ening is though to be due to the gas shocked by the high 
velocity outflow. High J fluxes measurements lead to esti- 
mate a shocked gas temperature between 1000 and 2000 
K, a H2 density between 10^ and 10^ cm~^ and a CO col- 
umn density between 10^^ and 10^^ cm~^. Our model both 
account of observed fluxes and lines profiles. This study 
show the necessity of high resolution observations of Orion 
BN/KL, as ISO-LWS observations not allow the spectral 
separation of the different gas components. Future obser- 
vations with HIFI, onboard the Far Infrared Space Tele- 
scope (FIRST) will allow the spectral separation of the 
PDR and the plateau component, unresolved with ISO, 
and characterise more precisely the Orion BN/KL star 
forming region. 
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In order to check the agreement of our model with ob- 
served lines profiles, we compared theoretical emission lines 
profiles of the three component with observed profiles. 
Theoretical profiles were obtained in convolving gaussian 
profiles by the PSF of the instrument. The line intensi- 
ties were inferred from lines fluxes predicted by our LVG 
model. The PSF was obtained by Vastel et al. in observ- 
ing thin [OI],[OIII] and [CII] lines in NGC7023 and G 
0.6-0.6. The theoretical profiles, superposed on observed 
lines, are shown on fig. [S] The theoretical profils are in 
good agreement with the observed ones. The predomi- 
nance of shocked gas emission at J^j, — 32 can explain 
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